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Abstract: The CEM2k and LAQGSM codes have been recently developed at Los Alamos National Laboratory to simulate 
nuclear reactions induced by particles and nuclei for a number of applications. We have benchmarked our codes against 
most available measured data at projectile energies from 10 MeV/A to 800 GeV/A and have compared our results with 
predictions of other current models used by the nuclear community. Here, we present a brief description of our codes and 
show illustrative results obtained with CEM2k and LAQGSM for A+p and A+A spallation, fission, and fragmentation 
reactions measured recently at GSI compared with predictions by other models. Further necessary work is outlined. 

Introduction 

During recent years, for a number of applications like Accelerator Transmutation of nuclear Waste (ATW), 
Accelerator Production of Tritium (APT), Rare Isotope Accelerator (RIA), Proton Radiography (Prad), astrophys- 
ical work for NASA, and other projects, we have developed at the Los Alamos National Laboratory an improved 
version of the Cascade-Exciton Model (GEM), contained in the code GEM2k, to describe nucleon-, pion-, and 
photon-induced reactions at incident energies up to about 5 GeV ^-[1] and the Los Alamos version of the Quark- 
Gluon String Model, realized in the high-energy code LAQGSM [7], to describe both particle- and nucleus- induced 
reactions at energies up to about 1 TeV/nucleon [7]-[Oj. 

Both codes have been tested against most of the available data and compared with predictions of other 
modern codes P-JSl- Our comparisons show that these codes describe a large variety of spallation, fission, and 
fragmentation reactions quite well and often have a better predictive power than some other available Monte-Carlo 
codes, thus they can be used as reliable event generators in different applications and in fundamental nuclear 
research. 

We have analyzed with CEM2k and LAQGSM all the A+p and A+A measurements done recently at GSI at 
energies near or below 1 GeV/nucleon for which we have results. The size of this paper allows us to present only 
a brief description of our models and a few results for proton-nucleus and nucleus- nucleus spallation, fission and 
fragmentation reactions measured at GSL Results for other reactions may be found in Refs. [T]-[T^. 

CEM2k and LAQGSM Codes 

A detailed description of the initial version of the CEM may be found in ^j, therefore we outline here only 
its basic assumptions. The CEM assumes that reactions occur in three stages. The first stage is the IntraNuclear 
Cascade (INC) in which primary particles can be re-scattered and produce secondary particles several times prior 
to absorption by or escape from the nucleus. The excited residual nucleus remaining after the cascade determines 
the particle-hole configuration that is the starting point for the preequilibrium stage of the reaction. The subse- 
quent relaxation of the nuclear excitation is treated in terms of an improved Modified Exciton Model (MEM) of 
preequilibrium decay followed by the equilibrium evaporative final stage of the reaction. Generally, all three stages 
contribute to experimentally measured outcomes. 

The improved cascade-exciton model in the code CEM2k differs from the older CEM95 version ^| by incor- 
porating new approximations for the elementary cross sections used in the cascade, more precise values for nuclear 
masses and pairing energies, a corrected systematics for the level-density parameters, adjusted cross sections for 
pion absorption on quasi-deuteron pairs inside a nucleus, the Pauli principle in the preequilibrium calculation, 
and an improved calculation of fission widths. Significant refinements and improvements in the algorithms used 
in many subroutines lead to a decrease of computing time by up to a factor of 6 for heavy nuclei, which is very 
important when performing simulations with transport codes. Essentially, GEM2k has a longer cascade stage, less 
preequilibrium emission, and a longer evaporation stage with a higher initial excitation energy, compared to its 
precursors GEM97 ilGj and CEM95 Besides the changes to CEM97 and CEM95 mentioned, we also made 

a number of other improvements and refinements, such as: (i) imposing momentum-energy conservation for each 
simulated event (the Monte-Carlo algorithm previously used in CEM provided momentum-energy conservation 
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only statistically, but not exactly for the cascade stage of each event), (ii) using real binding energies for nucleons 
at the cascade stage instead of the approximation of a constant separation energy of 7 MeV used in previous 
versions of the CEM, (iii) using reduced masses of particles in the calculation of their emission widths instead of 
using the approximation of no recoil used previously, and (iv) a better approximation of the total reaction cross 
sections. On the whole, this set of improvements leads to a much better description of particle spectra and yields 
of residual nuclei and a better agreement with available data for a variety of reactions. Details, examples, and 
further references may be found in ^-IH]- 

The Los Alamos version of the Quark-Gluon String Model (LAQGSM) ,7, is a further development of the 
Quark-Gluon String Model (QGSM) by Amelin, Gudima, and Toneev (see ^Tj and references therein) and is 
intended to describe both particle- and nucleus-induced reactions at energies up to about 1 TeV/nucleon. The core 
of the QGSM is built on a time-dependent version of the intranuclear-cascade model developed at Dubna, often 
referred in the literature simply as the Dubna intranuclear Cascade Model (DCM) (see ^Sj and references therein). 
The DCM models interactions of fast cascade particles ( "participants" ) with nucleon spectators of both the target 
and projectile nuclei and includes interactions of two participants (cascade particles) as well. It uses experimental 
cross sections (or those calculated by the Quark-Gluon String Model for energies above 4.5 GeV/nucleon) for these 
elementary interactions to simulate angular and energy distributions of cascade particles, also considering the Pauli 
exclusion principle. When the cascade stage of a reaction is completed, QGSM uses the coalescence model described 
in JH| to "create" high-energy d, t, '^He, and ''He by final-state interactions among emitted cascade nucleons outside 
of the colliding nuclei. After calculating the coalescence stage of a reaction, QGSM moves to the description of 
the last slow stages of the reaction, namely to preequilibrium decay and evaporation, with a possible competition 
of fission using the standard version of the CEM If the residual nuclei have atomic numbers with A < 13, 

QGSM uses the Fermi break-up model to calculate their further disintegration instead of using the preequilibrium 
and evaporation models. LAQGSM differs from QGSM by replacing the preequilibrium and evaporation parts of 
QGSM described according to the standard CEM with the new physics from CEM2k and has a number of 
improvements and refinements in the cascade, coalescence, and the Fermi break-up models (in the current version 
of LAQGSM, we use the Fermi break-up model only for A < 12). A detailed description of LAQGSM and further 
references may be found in and in our later publications 0] El E| ■ 

Originally, both CEM2k and LAQGSM were not able to describe fission reactions and production of light 
fragments heavier than ^He, as they had neither a high-energy-fission nor a fragmentation model. Recently, 
we addressed these problems by further improving our codes and by merging them with the Generalized 

Evaporation Model code GEM2 developed by Furihata [201 • 

Our current versions of CEM2k and LAQGSM were incorporated recently into the MARS and LAHET trans- 
port codes and are currently being incorporated into MCNPX. This will allow others to use our codes as event- 
generators in these transport codes to simulate reactions with targets of practically arbitrary geometry and nuclide 
composition. 

Illustrative Results 

The Generahzed Evaporation Model code GEM2 of Furihata [517 merged with both our CEM2k and LAQGSM 
takes into account evaporation of up to 66 types of particles and light fragments (from n to ^^Mg) from excited 
compound nuclei, while most other evaporation models used in the literature consider evaporation of only 6 types 
of particles, from n to "^He. It is interesting to see how important this is when analyzing the recent GSI A-l-p and 
A-l-A measurements. Fig. 1 gives us a quick and clear answer to this question: It is not important at all. We 
see that calculations by CEM2k-t-GEM2 taking into account up to 66 types of evaporated particles and fragments 
almost coincide with similar results calculated considering only 6 types of evaporated particles for all products 
measured recently at GSI |23 for the reaction 800 MeV/A ^^^Au-|-p. Similar results were obtained for other GSI 
measurements. We do see a big difference between the results of these "66" and "6" calculations, but only for 
products with A < 28 and Z < 12, and such light products have not been measured at GSI. Note that the "66" 
calculations require about 7-8 times more computing time than the "6" ones. This means that if we study only 
spallation and fission products and are not interested in light fragments, we can consider evaporation of only 6 
types of particles and save the computing time, getting results very close to the ones calculated with the more 
time consuming "66" option. But if we need to describe correctly all products from a reaction, including the light 
fragments, we will need to use the "66" option. All results presented below (and in previous publications) were 
obtained using the "66" option. 

Fig. 2 shows an example of a proton-induced reaction calculated by LAQGSM-f GEM2 measured at GSI 
in inverse kinematics, i.e., as 1 GeV/A ^°'^Pb-|-p [^H- For comparison, results obtained with the transport code 
LAIIET3 1221 using the Liege intranuclear cascade (INC) model INCL by Cugnon et al. 12^ merged with the GSI 
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Figure 1. The measured mass and charge 
distributions of the product yields from the 
reaction 800 MeV/A ^^^Au+p and of the 
mean kinetic energy of these products, and 
the mass distributions of the cross sections 
for the production of thirteen isotopes with 
the charge Z from 20 to 80 (circles) compared 
with CEM2k+GEM2 calculations taking into 
account evaporation up to 66 types of particles 
and light fragments (blue lines) and considering 
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evaporation/fission model ABLA by Schmidt 
et al. 25 , as incorporated into LAHET3 by 
J.-C. David of CEA, Saclay (see some details 
in Appendix 2). 

Comparison of our CEM2k+GEM2 and 
LAQGSM+GEM2 results with the third 
and the last proton-induced reaction mea- 
sured at GSI in inverse kinematics we found 
tabulated data, namely 1 GeV/A ^^^U-l-p 
[inillTl is shown in Figs. Al.l-Al.4 of Ap- 
pendix 1. For convenience, we divide the 
results for the products shown in these fig- 
ures into three groups: spallation, fission, 
and fragmentation. We compare our re- 
sults with calculation by LAHET3 using 
the INCL-hABLA |21I2^ option mentioned 
above, as well as using the ISABEL .23. in- 
tranuclear cascade model coupled with the 
Dresner evaporation code and the RAL 
fission model of Atchison (30] , as well as us- 
ing the Bertini INC |3J coupled with Dres- 
ner and Atchison [30] models. 
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Figure 2. Experimental [22] mass distributions of the cross 
sections of thirteen isotopes with the charge Z from 22 to 82 
compared with our LAQGSM-I-GEM2 calculations and results by 
LAHET3 |28| using the Cugnon et al. intranuclear cascade model 
INCL |24| merged with the GSI evaporation/fission model ABLA 
by Schmidt et al. j^. 

We performed all calculations of this reaction in 2002 (ex- 
cept with INCL-I-ABLA, see details in Appendix 2), after the 
measured spallation product cross sections were published in 
|26j . and published part of these results in a 2002 LANL Theo- 
retical Division Report of Activity "3^ . The experimental data 
on fission and fragmentation products were published only in 
2003 [13; therefore the CEM2k-|-GEM2, LAQGSM-hGEM2, 
IS ABEL+Dresner/ Atchison, and Bertini+Dresner/ Atchison 
results for fission and fragmentation products shown in Figs. 
Al.l-Al.4 are pure predictions, many of which agree amazingly 
well with the experimental data. 
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Figure 3. Comparison of preliminary experimental data taken by us 
from Fig. 5 of Carmen Villagrassa j33) on mass distribution of the yields 
of eight isotopes from Na to Mn produced in the reaction 1 GeV/A 
'^'^Fe + p (circles) with our CEM2k+GEM2 results and with calculations 
by LAHET3 using the INCL+ABLA and IS ABEL+Dresner/ Atchison 
options (lines), respectively. 

As do all Other models, CEM2k+GEM2 and LAQGSM+GEM2 
have many parameters. But all these parameters are fixed and all 
the results shown in the figures of this paper were calculated within 
a single approach, without fitting any parameters: We changed only 
the values of the mass and charge numbers of the projectile and target 
nuclei and the incident energy of the projectile in the input files of 
our codes. 

Fig. 3 shows the last proton-induced reaction we discuss in the 
present work, 1 GeV/A ^^Fe+p. In a way, our calculations shown in 
this figure can be considered also as predictions, as we do not have 
numerical values for the GSI measurements of this reaction and the 
circles shown in Fig. 3 as "GSI data" are taken by us from Fig. 5 of 
Carmen Villagrassa |33) and are only preliminary. 

In the rest of the paper we focus at nucleus-nucleus reactions mea- 
sured recently at GSI, and we start our analysis with the lightest 
target, d, namely with the reaction ^"^Pb(l GeV/A) -I- d [31] shown 
in Fig. 4. 
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Figure 4. Experimental |84| mass distributions of the cross sections 
of ten isotopes with the charge Z from 35 to 82 (circles) produced in 
the reaction 1 GeV/A ^"^Pb + d compared with our LAQGSM+GEM2 
results and with calculations by LAHET3 using the INCL+ABLA and 
IS ABEL-|-Dresner/ Atchison options (lines), respectively. 

One can see that LAQGSM-I-GEM2 describes quite well both the 
spallation and fission product cross sections and agrees with most of 
the GSI data with an accuracy of a factor of two or better similarly 
to the results of the INCL-f ABLA and IS ABEL+Dresner/ Atchison 
models. 

Fig. Al.5 in Appendix 1 shows another reaction on d: 1 GeV/A 
238y_j_j_ Only the yields of spallation products from this reac- 
tion measured at GSI are available to us In Fig. Al.5, 
we compare our LAQGSM+GEM2 results with all published data 
and with calculations by LAIIET3 using the INCL+ABLA and IS- 
ABEL+Drcsncr/ Atchison options. We see that as in the case of 
Pb+d, LAQGSM+GEM2 describes the U+d data quite well, bet- 
ter than do the the INCL+ABLA and IS ABEL+Dresner/ Atchison 
models. 

Our CEM2k+GEM2 and LAQGSM+GEM2 codes describe cor- 
rectly not only the cross sections of the spallation, fission, and frag- 
mentation products from various reactions, but also their mean ki- 
netic (recoil) energy. One example on this is shown in Fig. A1.6 of 
Appendix 1. 

Fig. 5 shows an example of a reaction on a heavier target, ^Be, 
namely the reaction 1 GeV/nucleon ^^Kr + ^Be measured by Voss 
[36', compared with our LAQGSM+GEM2 results. No fission mech- 
anism is involved in this reaction and all the measured products pub- 
lished in | 36| and shown in this figure are described by our code using 
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Figure 5. Comparison of all measured | 8(i| cross sections of products 
from the reaction *®Kr + ^Be at 1 GeV/nucleon (symbols) with our 
LAQGSM+GEM2 results (lines). 

only spallation. Although LAQGSM+GEM2 underestimates signif- 
icantly the yields of neutron-rich Rb isotopes, otherwise there is a 
good agreement between the calculations and data for all the other 
measured cross sections. 

Fig. 6 shows an example of a reaction on a heavier target, ^^Al, 
namely the reaction 790 MeV/nucleon ^^^Xe + ^^Al measured at 
GSI by Reinhold et al. [37] and compared with LAQGSM-I-GEM2 
results. Although both the projectile and target are heavier than 
for the example shown in Fig. 5, LAQGSM-t-GEM2 describes all the 
products from the reaction shown in Fig. 6 as well using only spalla- 
tion. A very good agreement between the data and calculations may 
be seen for all measured cross sections, except for the neutron-rich 
Cs isotopes, whose charge is bigger than that of initial Xe nuclei of 
the beam, being produced by picking up a proton from the Al target 
rather than by spallation processes. The situation observed in Fig. 
5 for the production of neutron-rich Rb isotopes involves the same 
process. 

Figure 6. Comparison of all measured |37| cross sections of products 
from the reaction ^^^Xe -1- ^^AI at 790 MeV/nucleon (filled circles) with 
our LAQGSM-I-GEM2 results (open circles). Isotopes from Fe to Y are 
not measured yet and we present here our predictions. 

Finally, Fig. Al.7 (Appendix 1) shows a heavy-ion-induced reaction 
measured at GSI |3Hll3Pj namely the yields of measured spallation 
products from the interaction of a 950 MeV/nucleon ^38^ beam with 
copper compared with our results. LAQGSM-I-GEM2 describes most 
of these data with an accuracy of a factor of two or better. 

Fig. A1.8 (Appendix 1) shows an example of several exotic re- 
actions, namely fragmentation of secondary beams of neutron-rich 
unstable 19^20,21q ^^^^ stable I'^'^^Q isotopes on ^^C targets at beam 
energies near 600 MeV/nucleon measured recently at GSI 001 j com- 
pared with our LAQGSM-I-GEM2 results. The secondary beams 
of 17-210 

ions were produced in the fragmentation of a primary 
^"Ar beam at 720 MeV/nucleon on a beryllium target (see more 
details in 001). A detailed discussion of our results on this reac- 
tion and more examples of nucleus-nucleus reactions analyzed with 
LAQGSM-HGEM2 may be found in [H]. 
From the results presented here and in the cited references, we conclude that CEM2k and LAQGSM describe well 
(and without any refitted parameters) a large variety of medium- and high-energy nuclear reactions induced both 
by nuclei and particles and are suitable for evaluations of nuclear data for applications and to study basic problems 
in nuclear reaction science. Merging our CEM2k and LAQGSM with the code GEM2 by Furihata allows us to 
describe reasonably well many fission and fragmentation reactions in addition to the spallation reactions already 
described well by CEM2k and LAQGSM. This does not mean that our codes are without problems. For instance, 
LAQGSM-I-GEM2 does not reproduce well the mass distributions for some fission-fragment elements from the 
reaction 1 GeV/A ^38^ _|_ 208p|-| measured recently at GSI |42j, although it still reproduces very well the integrated 
mass- and charge-distributions of all products. We think that the main reasons for this problem are the facts that 
the current version of LAQGSM does not take into account electromagnetic-induced fission '53' , and because the 
GEM2 code by Furihata merged at present with our LAQGSM does not consider at all the angular momentum 
of emitted particles, and of the compound nuclei. Both these factors are especially important for reactions with 
heavy ions and less important for reactions with light ions or protons; this would explain why the code works well 
in the case of reactions induced by particles and light and medium nuclei but fails in the case of U-|-Pb. Our work 
on CEM2k and LAQGSM is not completed; we continue their further development and improvement. Some details 
of our present work and plan for future may be found in Refs. [121 119| . 
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Figure A1.2. Comparison of measured |2E1 spallation product cross sections of the reaction ^^^U(l 
GeV/A) + p (filled circles) with results by LAQGSM+GEM2 and by LAHET3 using the INCL+ABLA, IS- 
ABEL+Dresner/ Atchison, and Bertini+Dresner/ Atchison options (open circles), respectively. Experimental data 
for isotopes from Tb to Ta and for Np and Pu are not yet available so we present here only our predictions. 
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Figure A1.3. Comparison of measured "ST fission-product cross sections of the reaction ^■^'*U(1 GeV/A) + p (filled 
circles) with results by LAQGSM+GEM2 and by LAHET3 using the INCL+ABLA, ISABEL+Dresner/Atchison, 
and Bertini -|-Dresner/ Atchison options (open circles), respectively. All our calculations (except for the 
INCL+ABLA, see details in Appendix 2) were done in 2002 and were published partially in a 2002 LANL Activity 
Report 121] while the data (27j were published and become available to us only in 2003. 
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Figure A1.4. Comparison of measured j2Z| fragmentation (and fission) product cross sections of the reaction 
238U(1 GeV/A) + p (filled circles) with results by LAQGSM+GEM2 and by LAHET3 using the INCL+ABLA, 
ISABEL+Dresner/ Atchison, and Bertini+Dresner/ Atchison options (open circles), respectively. Experimental data 
for isotopes from B to Co arc not yet available so we present here only our predictions. All our calculations (except 
for the INCL+ABLA, see details in Appendix 2) were done in 2002 and were published partially in a 2002 LANL 
Activity Report 32 while the data [21] were published and become available to us only in 2003. 
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Figure A1.5. Comparison of measured spallation-product cross sections of the reaction ^'^^U(l GeV/A) 
+ d (filled circles) with results by LAQGSM+GEM2 and by LAHET3 using the INCL+ABLA and IS- 
ABEL+Dresner/ Atchison options (open circles), respectively. Experimental data for Pu isotopes are not yet 
available so we present here only our predictions. 



11 



^ 110 

> 

^ 90 



A 70 r 



N 



50 r 



LU 

V 30 



-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 

1 A GeV '''U + p 




OGSI data 

— CEM2k+GEM2 

— LAQGSM+GEM2 



25 



35 45 55 65 

Charge number, Z 




125 



145 165 185 205 

Mass number, A 
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Figure A1.8. Cross sections of projectile fragments with nuclear charges Zf (shown on the top) and masses Af 
(shown on the bottom) produced from ^^~^^0 beams (shown on the right) in a ^^C target. Experimental data 
(filled circles) are from 00]. Open circles show our LAQGSM-I-GEM2 results for the measured cross sections and 
predictions for several unmeasured isotopes. Dashed lines show results by the the abrasion- ablation model 01] 
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Appendix 2 



We apologize for showing at the TRAMU workshop results calculated with the transport code LAHET3 [231 using 
the Liege intranuclear cascade model INCL by Cugnon et al. (24) and citing them as calculated using the evapora- 
tion/fission model ABLA by K.-H. Schmidt et al. |2S|, i.e., INCL+ABLA. The truth is that one of us performed 
these calculations with INCL using the default options of LAHET3 for evaporation and fission, which selects the 
Dresner evaporation model ^22, in conjunction with the Atchison fission model QQ^ but not ABLA. We thank Dr. 
Sylvie Leray to calling our attention to this confusion. To correct ourselves and to see how using different evapora- 
tion/fission models coupled with the same intranuclear cascade (INC) model affects the final results we performed 
after the TRAMU workshop calculations with LAIIET3 using the right INCL+ABLA option. In the following three 
pages, we present several figures with mass and charge distributions of the products from all reactions we discussed 
at TRAMU calculated with INCL+ABLA compared with similar results by INCL+Dresner/ Atchison showed at 
TRAMU (and cited mistakenly as INCL+ABLA) and by Bcrtini+Dresner/ Atchison, IS ABEL+Dresner/ Atchison, 
CEM2k+GEM2, and LAQGSM+GEM2 models discussed in our talk at TRAMU. One may see a big differ- 
ence between INCL+ABLA and INCL+Dresner/ Atchison results. On the whole, INCL+ABLA provides a better 
agreement with the measured fission products than when using the Dresner/ Atchison option to calculate evap- 
oration and fission, but strongly underestimates the spallation products, especially with masses near the bor- 
der between the spallation and fission regions, and agrees much worse with the data in this mass region than 
INCL+Dresner/ Atchison does for these products. We see that different evaporation/fission models coupled with 
the same INC model provide significantly different results in both the spallation and fission regions of products, 
suggesting that development of a reliable and universal evaporation/fission model is of a first priority for any 
transport code, independently of what INC model is used to describe the intranuclear-cascade stage of a reaction. 
Finally, we would like to warn future users of LAHET3 that when choosing the INCL option for the intranuclear 
cascade model in LAIIET3, the default options for evaporation will be Dresner and for fission, Atchison, and one 
needs to specify explicitely the option ABLA to get results by INCL+ABLA, otherwise LAIIET3 will provide 
results by INCL+Dresner/ Atchison, as happened to us before TRAMU. 
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Figure A2.1. Comparison of measured mass (left panel) and charge (right panel) distributions of the nuclides 
produced in the reactions 1 GcV/A ^"^Pb + p and 800 MeV/A ^^'^Au + p with results by LAHET3 using the 
Bertini+Dresner/ Atchison, IS ABEL+Dresner/ Atchison, INCL+Dresner/Atchison, and the INCL+ABLA options, 
and by the CEM2k+GEM2 and LAQGSM+GEM2 codes (color lines), respectively. Experimental data for ^ospb 
are from Ref. (22| and for ^^^Au, from Ref. |21| . 
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Figure A2.2. Comparison of measured mass (left panel) and charge (right panel) distributions of the nuclides 
produced in the reactions 1 GeV/A ^^^U + p and 1 GeV/A ^^Fe + p with results from LAHET3 using the 
Bertini+Dresner/ Atchison, IS ABEL+Dresner/ Atchison, INCL+Dresner/Atchison, and the INCL+ABLA options, 
and by the CEM2k+GEM2 and LAQGSM+GEM2 codes (color lines), respectively. Experimental data for ^asy 
are from Refs. I26L I27| . The 1 GeV/A GSI data for ^^Fe were obtained using the systematics by Morrissey 0^ 
and the Goldhaber model 0^ and are extracted by us from Figs. 7 and 8 of Carmen Villagrassa earlier Fe 
data at 1.05 GeV/A by Zeitlin et al. gSI, 1.086 GeV/A by Webber et al. and 1.88 GeV/A by WestfaU et al. 
0S| (all tabulated in Tab. 7 of Ref. 0^) are also shown for comparison. 
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Figure A2.3. Comparison of measured mass (left panel) and charge (right panel) distributions of the nuclides 
produced in the reactions 1 GeV/A ^^^Pb + d and 1 GeV/A ^^^U + d with results from LAHET3 using the IS- 
ABEL+Dresner/ Atchison, INCL+Dresner/Atchison, and the INCL+ABLA options, and by our LAQGSM+GEM2 
(color lines), respectively. Experimental data for ^"*Pb are from Tabs. 2 and 3 and Fig. 13 of Ref. Pl] and for 
from Ref. 
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